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We briefly review the recent experimental results on possible variations of the fine structure con- 
\j^ • stant a on the cosmological time scale and its position dependence. We outline the theoretical 

^P . grounds for the assumption that a might be variable, mention some phenomenological models 

incorporating a variable a into the context of modern cosmology and discuss the significance of 

possible a variations for theoretical and practical metrology. 
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1. Introduction 

The idea that fundamental physical constants (FPC), including the Planck constant, the 
Q I speed of light, stable particle masses and the coupling constants of the basic interactions 

O ' (above all, Newton's gravitational constant G and the fine structure constant a) may be 

Q ■ variable, was put forward for the first time by Dirac in 1937, soon after the advent of Fried- 

^ . mann's expanding Universe models (1922) and their confirmation by Hubble (1929). In the 

--«^ ! subsequent decades, the possible FPC variability in the course of cosmological evolution as 

qh| well as their dependence on position in space, on the magnitude of physical fields or on 

spatial and energy scales occupied quite a prominent place in theoretical physics and cos- 
W)' mology, including many works of the VNIIMS theoretical group, Refs. [1-11] and others. In 

addition to the above-mentioned FPC, there appeared other "candidate variables" such as, 
K> ! for instance, the cosmological constant and the coupling constants of the weak and strong 

H I interactions. The basic motivation for such studies has been, and remains to be, the idea of 

unity of physical interactions. If the conjectured unified interaction splits in different physical 
conditions into specific interactions known from the experiment, the FPCs that character- 
ize these interactions should depend on the circumstances of their manifestation, above all, 
on the spatial and energy scales and on the fundamental field intensities (a recent detailed 
discussion may be found in Ref. [10]). 

Thus, it is quite common to believe that unification of all physical interactions requires 
a number of additional (so-called internal) space-time dimensions. At present, there is no 
well established way of transition from multidimensional physics to the conventional four- 
dimensional picture. It is, however, clear that such a transition requires some procedures 
(such as dimensional reduction and compactification) which connect the four-dimensional 
effective constants with true multidimensional constants and the characteristics of extra di- 
mensions. Thus, their size may suffer a non-trivial evolution in time and may depend on the 
position in external space. 

As to observational data on FPC variation, there have been very few positive results till 
the recent years. Thus, some researchers pointed toward a possible G variability according 
to astronomical observations, but, to our knowledge, neither of these results were confirmed 
afterwards (see Ref. [10] for detail). However, the very end of the 20th century was marked 



with the advent of observational evidence [12] on the variabihty of the electromagnetic cou- 
pling constant, the fine structure constant 

a = e^ /{he) (Gaussian units), or 

a = e^/(47r5o^c) (SI units), (1) 

where e is the electron charge, c is the speed of light in vacuum and Eq is the permittivity of 
free space. This and subsequent works have made the subject of FPC variability, and above 
all the variability of a, highly topical. This subject is no longer regarded purely theoretical, 
and nowadays it rapidly increases its experimental basis. 

The electromagnetic interaction, whose intensity is characterized by the constant a, is of 
primary importance for the macroscopic structure of matter and in an overwhelming majority 
of the observed phenomena. It is for this reason that a possible variability of a would lead to 
variability of the majority of existing references of physical quantities, including such basic 
references as those of length, mass, time/frequency. 

The constant a, whose CODATA-recommended value is 

a = 7.297362533(27) x lO^^ ^ 1/137.03599 

(the figures in parentheses mean the uncertainty in the last two meaningful digits), is one of 
the most precisely measured FPC: the relative uncertainty is about 3.7 x 10^^ . Moreover, its 
dimensionless nature makes it independent of any conventions or systems of physical units. 
In what follows, we will briefly describe the basic experimental and observational data on 
the possible a variability and outline the main characteristics of theoretical models which, 
to one or other extent, conform to the observations. In the last section, we formulate some 
conclusions and discuss the signiflcance of this research area for theoretical and practical 
metrology. 

2. Empirical data on a variability 

A great number of relevant observations and experiments are described and discussed in the 
comprehensive review [13] and in the shorter articles [14,15]. In this section, we mention 
only some of the data presented there, seeming to be the most convincing, and add some 
estimates and inferences of the last two years. 

2.1. Laboratory data 

The laboratory bounds on a variability are obtained by comparing the long-term behaviour of 
oscillators with frequencies depending on a in different ways. In practice, such a comparison 
leads to experimental bounds on certain combinations of FPCs. 

To our knowledge, the flrst comparisons of this kind were performed in 1974-76 by two 
groups independently. Turneaure and Stein [16], at Stanford University, used superconduct- 
ing microwave cavities with a quality factor of about 4 x 10^'^ and a resonant frequency ~ 8.6 
GHz. The frequencies of cavity-controlled oscillators were compared over a period of 12 days 
with a caesium-beam atomic clock, and the relative drift was (—0.4 ±3.4) x 10"^^ /day. This 
led to the conclusion [16] that a possible annual variation of the product 

gp ■ {rrie/M) ■ a 



was less than 4.1 x 10^^^ , where Qp is the proton gyromagnetic ratio, nie is the electron 
mass and M is the Cs nuclear mass. 

Kolosnitsyn et al. [11], at VNIIFTRI, the Russian time-frequency standard keeper, com- 
pared the behaviour of molecular clocks operating with ammonia molecular beams and a 
group of caesium frequency standards. The caesium standard formed a scale for measuring 
the frequency of the ammonia generators. The results of measurements for six years (1964- 
1969) led to a number of estimates for a/a, all of them constraining this variation by a few 
units times 10^"*^^; the best estimate was 

\a/a\ < 1.1 X lO^^yr^^ 

Other early results on the stability of a and e are described in Petley's book [17]. 

The recent estimates of possible a variations have diminished the uncertainties by 4 
to 5 orders of magnitude as compared with [11] due to a considerably increased stability of 
frequency standards. Thus, according to [18], a comparison of hyperfine transition frequencies 
in ^'''Rb and ^^^Cs over a period of about 4 years has shown that 

dln{uj^^/ucs)/dt = (0.2 ± 7.0) x lO^^V^"^ 

at la. Neglecting possible changes in the the weak and strong interaction constants and thus 
in the nuclear magnetic moments, it leads to 

a/a = {-0 A ± 16) xlO-^^yr~\ (2) 

where a = da/dt. Another experiment [19] compared an electric quadrupole transition in 
i99jjg+ ^Q ^Yie ground-state hyperfine splitting of ^^^Cs over 3 years, showing that 

|dln(z/Hg/z/cs)/c^t| < 7.0 X lO^^^r^i 

and leading to 

|d/a| < 1.2 X lO-^^r-i (3) 

under the assumption that the gyromagnetic factor qq^ and rric/mp are invariable. 

Ref. [20] reported a comparison of the absolute IS* — 25* transition in atomic hydrogen 
to the ground state of caesium, which, combined with the results of Refs. [11, 19], yielded the 
constraint 

d/a = (-0.9±2.9) X lO^^Vr"^- (4) 

A comparison [21] of optical transitions in ^^^ Yb+ with a caesium standard for 2.8 years 
has shown that dln{i'Y\^/i'cs)/dt = (—1.2 ± 4.4) x 10~^^yr~^ which leads to 

d/a = (-0.3±2.0) X lO^^^r-^ (5) 

The laboratory measurements thus lead to results which, as follows from Ref. [20], may 
be combined to yield very hard constraints for separate basic constants including a. The 
above estimates show that the allowed variations of a in the modern epoch are bounded to 
a few units of 10~^^ per year. 



Further progress in the near future is expected on the basis of improved frequency stan- 
dards as well as development of new methods of atomic experiments. One such method [22] 
is based on crossing of atomic levels in two-electron highly charged ions of different atomic 
numbers Z . It is claimed that the effect of possible variation of a becomes strongly (by a 
factor of about 1000) enhanced when studied near such crossing points. Another method [23] 
is to use atomic transitions between narrow lines that have an enhanced sensitivity to a pos- 
sible variation of a . The authors present a number of such transitions and claim that their 
method effectively suppresses the systematic sources of uncertainty that are unavoidable in 
conventional high-resolution spectroscopic measurements. 

An improvement of at least an order of magnitude in current a estimates can be achieved 
in space experiments planned for 2008-2009 within the PHARAO/ACES Project of the Euro- 
pean Space Agency [24,25]. PHARAO, the "Projet d'Horloge Atomique par Refroidissement 
d'Atomes en Orbite", which has started in 1993, combines laser cooling techniques and a 
microgravity environment in a satellite orbit, allowing the development of space clocks with 
unprecedented performance. Its declared objectives are, among others [24], (i) to achieve 
time and frequency transfer with stability better than 10~^^ and (ii) to perform fundamental 
physics tests. Its detailed description can be found in [25]. 

One can also mention quite a different type of a variability predicted by the electroweak 
theory which concerns strengthening of the electromagnetic coupling with increasing momen- 
tum transfer Q^ in particle interactions at energies approaching the electroweak unification 
energy [26,27]. This effect has been confirmed experimentally [26] by analyzing the processes 
e+e^ -^ /i+yU~ and e^e" -^ e~^e~fi^fi~ at the TRISTAN e+e^ collider at KEK (Japan). 
The data were accumulated at an average center-of-mass energy of 57.77 GeV, and a was 
measured to change from its known value of a~^ ~ 137.0 at Q^ = to 

a~^ = 128.5 ± 1.8 (stat) ± 0.7 (syst) at Q^ = (57.77 GeV/cf. 

The latter variation is physically interpreted as a result of quantum vacuum properties (see, 
e.g., [27] and references therein) and manifests itself at energies much higher than those 
relevant to atomic spectra. It has therefore nothing to do with temporal and spatial variations 
of a at astrophysical and cosmological scales. 

2.2. Geochemical data 

Very hard constraints on variations of a have been obtained from studies of the so-called 
Oklo phenomenon combined with data on the lifetimes of long-lived radioactive isotopes. 

The Oklo phenomenon is a natural nuclear reactor that operated for about 200,000 years 
approximately two billion years ago, which, regarding cosmological observations, corresponds 
to redshifts around z ~ 0.14. The products of its operation were discovered in 1972 at 
Oklo uranium mine in Gabon, West Africa. The isotopic abundances of the yields bears 
information on the nuclear reaction rates at that time, which in turn depended on the current 
values of a. One of the key quantities measured is the ratio g2^Sm/g|^Sm of two light 
isotopes of samarium which are not fission products. This ratio is of the order of 0.9 in 
normal samarium, but is about 0.02 in Oklo ores. Such a low value is explained by depletion 
of 62^Sm by thermal neutrons to which it was exposed while the reactor was active. The 



capture cross section of a thermal neutron by ^^^Sm, i.e., the reaction 

i^^Sm + n ^1^° Sm + 7 

has a resonant energy E^ — 0.0973 eV, which is a consequence of a nearly cancellation 
between electromagnetic and strong interaction forces [28]. A detailed analysis has shown [29] 
that, assuming that E^ varied only due to an a dependence of the electromagnetic energy, 
the following constraint is valid: 

Aa/a = (0.15 ±1.05) X 10"^ (6) 

at 2o" level. The accuracy of the method is explained by the tiny value of the resonant energy 
-Er ~ 0.1 eV compared to its sensitivity to a variation of a, dEr/dhia ~ —1 MeV, so that 
the expected variations are smaller than 0.1 eV/1 MeV ~ 10^^. 

It was later noticed [30] that two ranges of variations could be compatible with the Oklo 
data: 

Aa/a = (-0.8±1.0) X 10~^ Aa/« = (8.8 ± 0.7) x 10"^ (7) 

but the second range is hardly compatible with the isotopic ratio of gadolinium. Note that 
the first range, assuming that the time variations of a (if any) occur uniformly, is translated 
to the variation rate 

d/a = (-0.4 + 0.5)xl0-^^yr-\ (8) 

which is apparently the tightest of the existing constraints. 

Some estimates have been recently obtained under the additional assumption that the 
low energy neutron spectrum is well described by a Maxwell-Boltzmann distribution [31- 
33]. Taking into account the possible variation of the strange quark mass rris along with 
the assumption that all fundamental couplings vary independently led [33] to the allowed 
variation Aa/a < (1 — 5) x 10~^°, even more stringent than above. 

Another result [34], free of poorly grounded assumptions, seems more realistic: 

Aa/a = (4.5 + 1) X 10-^ d/a = (2.25 + 0.5) x 10~^^ yr-^ (9) 

The most recent a estimates from the Oklo phenomenon [35] have been made on the basis 
of modern methods of reactor physics, a detailed computer model of zone RZ2 of the Oklo 
reactor and full-scale calculations using two different Monte Carlo codes: the Russian code 
MCU REA and the worldwide known code MCNP (USA). Both codes have produced close 
results. It was claimed that many details of these calculations differed from the previous 
work (e.g., the averaged cross-section of Sm and its dependence on the shift of resonance 
position due to variation of fundamental constants). Still no evidence on a variations was 
found, with the resulting constraint 

-4 X lO^^V"^ < {da/dt)/a < 3 x lO^^V"^- 

A further improvement in the accuracy of these limits is promised. 

Some conclusions have also been made by analyzing the readioactive decay of nuclides 
with known a dependence of the decay rate. The best constraint was obtained from beta 



decay of rhenium into osmium with electron emission: as noted by Peebles and Dicke [36] , 
an extremely low value of the decay energy, about 2.5 keV, makes it very sensitive to a 
variations. Still the constraints obtained are much weaker than (jH)) or ©• 
Laboratory studies of meteorites have led to the result [37] 

Aa/a = (8±16)x 10-^ (10) 

for the recent 4.5 billion years, which corresponds to z ~ 0.45. Translated to the changing 
rate (again assuming its uniformity) gives approximately 

d/a = (2 ± 4) X lO-^Vyr- (H) 

A shortcoming of this method is its indirect nature: the Re/Os ratio is measured in iron 
meteorites whose age is not determined directly. Besides, the constraint (fTIH) only concerns 
the mean value of a for 4.5 billion years rather than its instant variation rate. 

2.3. Astrophysical data 

The first and not very confident observational indication on a possible variability of a has 
appeared from studies of remote quasars. The observed quasar absorption spectra were 
compared with the corresponding laboratory spectra. The details of the method are described 
in Refs. [12,38-41] and others. 

One of the recent results [39] has been obtained by analyzing the lines of SilV in 15 
systems and improves the previous estimates by a factor of three: 

Aa/a = (0.15 ±0.43) X 10-^ 1.59 < z < 2.92. (12) 

This result does not confirm a variations. 

The many-multiplet method of Webb et al. [12,40] rests on a comparison of a few lines 
in different samples of sources: one compares the shifts of lines which are sensitive and 
insensitive to a variations. A recent result [20] obtained from a study of 128 systems in the 
range 0.5 < z < 3 points at lower values of a in the past: 

Aa/a = (-0.54 ±0.12) X 10-^ 0.5 < z < 3, (13) 

thus confirming the previous conclusions [12,40]. The authors did not find systematic errors. 
There are, however, results directly opposite to these. Thus, according to [42], 

Aa/a = (-0.1 ± 1.7) X 10-^ ^ = 1.15 (14) 

(from an analysis of absorption lines of iron of a single quasar). Similar results have been 
found for a system of absorption lines of the quasar Q 1101-264 [43] with the redshift z = 
1.839: 

Aa/a = (2.4±3.8) X 10~^ ^ = 1.839. (15) 

A unified sample of Fe II lines gave Aa/a = (—0.4 ± 1.5) x 10~^ for two systems with 
z = 1.15 and 1.839. 



The use of radio and millimeter quasar absorption spectra gave an independent estimate 
(unfortunately, for comparatively small redshifts only) 

|Aa/a| < 8.5 X 10-^ z = 0.25 - 0.68. (16) 

As noted by Barrow [44], the same quasar absorption spectra analysis leads to an upper 
bound on spatial variations of a like |Aa|/a < 3 x 10~^ at 3 Gpc distance scale because of 
a wide distribution of the target absorption systems over the sky. 

An alternative method is to study emission rather than absorption spectra. There are very 
few such estimates since this method, being sufficiently simple and direct, is less sensitive and 
is harder to apply to sources with large redshifts. Recent measurements of strong emission 
lines of O III in a sample of 165 quasars have led to the result [45] 

Aa/a = (1.2±0.7) xlO-^ z = 0.16 - 0.8. (17) 

Essential improvements in emission line measurements and analysis is anticipated in the near 
future [46], partly owing to the use of many-multiplet methods already employed in the 
absorption spectra analysis. 

The estimate ()17|) as well as tentative estimates of the Cambridge group (England) [47] 
point at higher values of a in the past whereas the results related to absorption lines give 
preference to smaller a in the past. It appears tempting to conclude that this apparent 
discrepancy may be explained by some unidentified systematic errors. Arguments in favour 
of this viewpoint have been put forward by Bandeira and Corbelli [48], who found that at 
least the absorption line analysis could be subject to systematic effects, related to estimation 
of different sets of atomic transitions at different redshifts, and these hidden correlations thus 
could mimic a variable a. 

Another, "optimistic" viewpoint is, however, possible [47], that this is simply an indication 
of spatial variations of a since the emission and absorption methods are sensitive to the values 
of a in drastically different environments. 

One can conclude that the present astrophysical data on a variation are rather contra- 
dictory and need further verification and improvement. 

2.4. Cosmological data 

The observed anisotropy of the cosmic microwave background (CMB) and the abundance of 
light elements formed in primordial nucleosynthesis (PNS) constrain the variations of a on 
the cosmological scale on the level 10~^ . These constraints concern much greater redshifts z 
than any others: z ~ 1000 for CMB and z ~ 10^° for PNS. 

Variations of a change the Thomson scattering cross-section and hence the transparency 
of the medium and finally the recombination time. 

A recent analysis of the WMAP (Wilkinson Microwave Anisotropy Probe) data on CMB 
anisotropy has led to the result [49] 

Aa/a = (-1.5 ±3.5) X 10-2, (18) 

for z ~ 10^. It is, however, well known that a similar effect on CMB anisotropy could be due 
to variations of the gravitational constant [50]. Besides, consistent constraints on a variations 
of the order of 1% could only be obtained from CMB analysis combined with independent 



data on the cosmological parameters characterizing the expansion of the Universe at times 
close to the recombination epoch. 

Barrow [44] used the observed CMB isotropy at large angular scales to derive strong 
observational limits on any possible large-scale spatial variation in the values of a and other 
FPC (the electron to proton mass ratio and the Newtonian gravitational constant) assuming 
that their space-time evolution is driven by a scalar field. The constraints are strongly 
model-dependent. Thus, large-scale spatial fluctuations of a are bounded by 2 x 10~^ in the 
BSBM theory (see below) and by 1.2 x 10^*^ in the varying speed of light theories. These 
derived bounds are significantly stronger than any others, obtainable by direct observations 
of astrophysical objects. It should be noted, however, that these bounds concern variations 
between mutually remote regions of space with approximately equal physical conditions and 
do not apply to the possible discrepancy between the results of emission and absorption 
spectra of quasars (see the end of the previous section): in the latter case, the emitters and 
absorbers are apparently characterized by drastically different gravitational fields, matter 
densities and temperatures. 

The PNS theory predicts light element formation in the early Universe, and their resulting 
abundances depend on a delicate balance between the Universe expansion and the weak 
interaction strength which in turn determines the proton to neutron number ratio at the 
beginning of PNS. Ultimately, the predicted light element abundances depend on a number 
of fundamental constants, a being only one of them. 

A recent study incorporating seven parameters [51] has yielded 

Aa/a = (6 ± 4) X 10-^ 2 ~ 10^°. (19) 

A considerable improvement of the above estimates is expected with the appearance and 
analysis of new observational data from the WMAP and Planck Surveyor satellites. 

According to Ref. [47], the CMB and PNS analysis on the whole favours slightly smaller 
values of a in the past compared to its modern value. The absence of variations also remains 
admissible on the level of la . 

The bounds (fTHj) and ()19|) cannot be directly translated into estimates of the admissible 
current values of a since, to do so, one has to assume a certain dependence a(t) which is not 
only different in different theories, but even in specific models of a given theory; see examples 
in Sec. 3. 

2.5. Comparing data of different origin 

The above discussion of possible FPC variability combined observational evidence from 
quasar spectra with Solar system data and laboratory constraints. Most of the studies im- 
plicitly assume that local and cosmological observations are directly comparable. This is, 
however, a strong assumption which may prove to be wrong. If a given FPC, say, a, depends 
on a scalar field (^, a slowly varying cosmological background of may substantially differ 
from its local values in a specific galaxy, stellar cluster or another gravitationally bound 
object. One can freely imagine, for instance, that the Oklo data only reflect the value of a 
stationary galactic scalar field whereas the quasar data feel cosmological effects. 

Shaw and Barrow [52] have found that it is probably not the case, and local variations of 
the "constants" are able to track their global variations. A construction involving matched 



asymptotic expansions within a sufficiently wide scalar field model is analyzed assuming 
that the scalar (dilaton) field is only weakly coupled to gravity and has a negligible effect 
on the background space-time geometry; a local inhomogeneity was characterized by the 
Tolman-Bondi class of spherically symmetric metrics describing the evolution of a dust cloud 
in the presence of a cosmo logical constant. It was concluded that local virialization does not 
stabilize the value of the dilaton and protect it from any global cosmological variation. This 
indicates that local terrestrial and Solar-system experiments really measure the effects of 
varying "constants" of Nature occurring on cosmological scales to computable precision [52] . 

Somewhat different results have been obtained by Mota and Barrow [53,54]. They studied 
the space and time evolution of a using the BSBM varying-a theory (see Sec. 3.1) and a 
spherical collapse model for cosmological structure formation, considered as a perturbation 
to a homogeneous and isotropic cosmological background. Different assumptions were used 
on the dark energy equation of state and on the coupling of a to the matter fields. It was 
found that, independently of the model of structure formation one considers, there is always 
a difference between the values of alpha in a virialized overdensity and in the background. In 
some models, especially at low redshifts, the difference depends on the time when virialization 
has occurred and the equation of state of the dark energy. At low redshifts, when the dark 
energy starts to dominate the cosmological expansion, the difference between alpha in a 
cluster and in the background grows. 

According to these studies, it may happen that a and a vary from galaxy to galaxy even 
at the same redshift, but the evolution of a is the same in different sub-galactic stellar systems 
belonging to a single galaxy, or even in galaxies belonging to the same galaxy cluster [55]. It 
is, however, clear that the problem should be studied more thoroughly in a wider range of 
models. 

More generally, one can conclude that the whole set of modern experimental and ob- 
servational data leads to very stringent limits of putative a variations but leaves open the 
question of their real existence. Even more uncertainty follows from the possible difference 
between local and global a variations and its model dependence. 

3. Theoretical models that predict varying a 

Despite the weak experimental status of varying a, there have been a great number of 
theoretical studies in the recent years, treating a as a function of certain physical fields 
which, in general, change from point to point in four- dimensional space-time. 

It should be noted that the researchers themselves stress a tentative nature of such con- 
structions. Thus, according to [47], "any model compatible with all empirical data at a 
given time, is certainly wrong since at any time some of the current data are erroneous". It 
is therefore useless to build artificial models with a large number of free parameters: it is 
highly probable that such a model will be incompatible with the newest data as soon as its 
description is published. 

In studies dealing with multidimensional theories, such as Kaluza-Klein theories and 
different versions of string theory (regarded at present as the most promising unification 
theories), the dimensionless constants of four- dimensional physics are shown to be dynamical 
quantities [56,57]. However, the specific form of such dynamics not only depends on the 
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version of the Kaluza-Klein or string theory, but also on the compactification scheme and 
the form of dilatonic couphng adopted within the chosen version of the theory. 

Accordingly, at the present stage, comparatively simple and natural phenomenological 
models are preferable, which, on the one hand, should conform to the idea of unity of all 
physical interactions and, on the other, should not only be able to predict the variability of a 
but also suggest solutions to other problems of modern cosmology, such as the cosmological 
constant problem and the dark matter and dark energy problems. We will briefly outline two 
such models and mention some other approaches. 

3.1. Bekenstein-Sandvik-Barrow-Magueijo (BSBM) theories 

This class of theories [58, 59] assumes constant values of the speed of light and Planck's 
constant; accordingly, a variability of a is achieved due to variability of the electron charge e 
or the permittivity of vacuum. Thus instead of e = const ane assumes e = e^eix) , where cq = 
const while e is a dimensionless scalar field depending on the space-time coordinates. Such 
theories preserve the usual properties of local gauge and Lorentz invariance and causality. 
In these theories, the conventional expression for the electromagnetic field tensor in terms of 
the 4-potential A^ is replaced by 

F^, = (1/5)[M,),^-M^),,], (20) 

which takes its usual form if e = const . The electromagnetic field action has its usual form 
while the scalar field dynamics is described by the usual kinetic term in the Lagrangian (writ- 
ten in the form e^^s'^ /e^). The theories of this class differ from one another in assumptions 
on how the e field is coupled to the rest of matter. 

The cosmological models obtained in these theories predict time variations of a (e.g., 
~ In t in the matter-dominated phase in a spatially flat cosmology) and spatial variations of 
a whose nature and evolution depends on the dark energy equation of state and the coupling 
of a to the matter fields [53,54]. Inclusion of the effects of inhomogeneity, so that the 
dependence a{t) differs in the cosmological background and in virialized overdensities and 
may naturally lead to no observable local time variations of a on Earth and in our Galaxy 
even though such variations can be significant on quasar scales. 

3.2. Fujii's scalar-tensor theory 

The theory suggested by Fujii [60] contains, in addition to the metric tensor, two scalar fields. 
One of them (a) is a dilatonic Brans-Dicke type field with the potential V{a) = Ae~^^'^ 
(in the Einstein conformal frame, or picture), where A is a constant playing the role of a 
cosmological constant in the Jordan picture; ( is the so-called non-minimal coupling constant 
characterizing a coupling between the a field and the metric in the Jordan picture. Another 
scalar field has a conventional kinetic term and interacts with a in such a way that their 
common potential in Einstein's picture has the form 

y(0, a) = e-^'^'^[A + lm^(j)\l + 7 sm{ka))], (21) 

where m, 7, k are constants (approximately of order unity in the Planck system of units). 

Considering Jordan's picture as a physical one, the author studies the modification of 
standard quantum electrodynamics in this theory caused by the scalar fields and concludes 
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that the effective electric charges are a-dependent. The time dependence of a itself is affected 
by the field which leads to decaying oscillations of a in the modern epoch. 

A cosmological model built in the framework of this theory predicts comparatively rapid 
changes of a in the epoch corresponding to 2; > 1 (but close enough to 1) and much slower 
changes in the modern epoch. This allows one to reconcile the strong bounds related to the 
Oklo phenomenon and the conclusions [12,40] on rather a rapid growth of a in the previous 
epoch. The field a is also able to play the role of the so-called quintessence explaining the 
Universe acceleration observed in modern cosmology. 

Other predictions of Fujii's theory are the existence of a non-Newtonian short-range com- 
ponent of gravity and a weak equivalence principle violation. The orders of magnitude of 
both effects are compatible with modern experimental constraints. 

3.3. Other models 

We have briefly described two sufficiently well-known and elaborated phenomenological the- 
ories. Let us also mention some alternative approaches. 

Anchordoqui and Goldberg [61] discussed the so-called quintessence, i.e., a minimally 
coupled scalar field playing the role of an "accelerator" of the Universe expansion, as a possible 
generator of a variability. They considered scalar field potentials combining a number of 
exponential functions with different coefficients. Like Fujii's theory, this model makes it 
possible to reconcile the tough constraints obtained from the Oklo reactor with the results 
of Webb et al. obtained from quasar absorption spectra: they become compatible due to 
slowing down of the scalar field evolution in approaching the present epoch. 

Damour and Polyakov [62] suggest a more general approach: among the models of the 
Universe containing ordinary matter and a real scalar field, one singles out a class of low- 
energy models possessing invariance under diffeomorphisms and gauge transformations. The 
resulting theories have a Born-Infeld type Lagrangian structure. In such theories, the masses 
and (yf-factors of fermions, as well as their electromagnetic coupling constants, are scalar field 
dependent, whereas their electric charges and gravitational coupling may remain constant. 
One then analyzes the effect of all variable factors on the observed atomic spectra. The 
latter, as follows from Ref. [62], bear information on the possible variability of a as well as 
the gf -factors and particle mass ratios, and, moreover, all predictions concerning the atomic 
spectra turn out to be conformal gauge independent. 

Bertolami et al. [63], among others, stress the effective Lorentz and CPT symmetry vio- 
lation that accompanies variable fundamental constants in any theory where these constants 
are scalar field dependent since variable scalar fields inevitably create preferred space-time di- 
rections. As en example of such a theory, they analyze N=4 supergravity in four dimensions, 
containing two scalar fields, and show that some of its solutions lead to a{t) dependence com- 
patible with the astrophysical data, though at the expense of some fine tuning. In another 
paper [64], Bento, Bertolami and Santos considered a cosmological model with a two-field 
quintessence, i.e., a set of two scalar fields and ■0, minimally coupled to the Einstein gravity 
and possessing the potential 

y(0,^)=e-^'^P(0,^), (22) 

where P(0, ip) is a third-order polynomial depending on both fields, while the electromagnetic 
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field Lagrangian had the form 

L,_^ = \B{<P,^)F>^''F^,, (23) 

B{(f),ilj) being a linear function of both scalars. It was found that some solutions for this 
system describe a cosmology with a transient late period of accelerated expansion, making it 
possible to fit the data arising from quasar absorption spectra and comply with the bounds 
on the variation of a from the Oklo reactor, meteorite analysis, atomic clock measurements, 
cosmic microwave background radiation, and big bang nucleosynthesis. 

Kirillov (see [65] and references therein) suggested a modified field theory which is aimed at 
explaining all effects ascribed to dark matter and dark energy without explicitly introducing 
them. A phenomenological manifestation of this theory may be expressed as a modification 
of the law of gravity, such that, for a point source of mass Mq, the gravitational potential 
reads 

= -^^[l + /(r)], (24) 



where the function /(r) may be chosen in the form 

f(r) = I *^^/^°'* ln(rinaxe/r), r < r^ax, ^25) 

the parameter tq, having the order of 1-5 kpc, is different for different galaxies while rmax 
only depends on cosmological time and is now of the order of 100 Mpc. At r > rmax, Newton's 
law is restored but with an effective mass M^ax = ^o(l + ^max/''"o)- It is shown that such 
a modification is able to account for all dark matter effects observed and, in particular, to 
explain the whole variety of "dark matter halos" in astrophysical systems. It is also shown 
that, in the modified theory, the galaxy formation process leads to an effective equaton of 
state of the fictituous "dark matter" with w < —1/3 {w being the ratio of pressure to energy 
density) and therefore it can play the role of dark energy. The standard picture of the early 
Universe is preserved, but the theory predicts variation of all interaction constants (including 
G and a) with the same time dependence [66]. 

4. Concluding remarks. Physical and metrological significance of 
FPC variations 

General theoretical considerations related to the necessary unification of all interactions and 
the requirement of unity of the physical science lead, probably inevitably, to the idea of a 
dynamical nature of all known FPCs or at least a greater part of them. Thus more surprising 
is that so much effort is needed to discover their variability than the belief that they are 
variable. Meanwhile, the experimental confirmations of these ideas are at present very poor 
and unreliable. This certainly concerns all FPCs, not only a. 

An important and natural feature of all models is the interrelation of the predicted varia- 
tions of different FPCs. Therefore, in both planning and interpretation of future experiments, 
such interrelations should be taken into account. Thus, in the context of string theories, it 
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turns out [62] that the quantum chromodynamics constant Aqcd and the weak interaction 
constant v are "even more variable" than a [15]: 

AAqcd/Aqcd ~ 30Aa/a, ^v/v ~ 80A«/«. (26) 

According to [67], any theoretical model, predicting cosmological variations of a over 
10^^ due to changes of some long-range scalar field, also predicts violations of the weak 
equivalence principle (WEP) over 10~^^. To discover such variations, it is only necessary 
to gain one order of magnitude as compared to the modern experimental constraints. This 
endows special significance to the space experiments for testing the WEP planned for the 
nearest years: MICROSCOPE (planned sensitivity 10"^^), STEP (up to 10"^^ [68]) and 
others. 

Thus new space experiments and observations in the area of gravitation, cosmology and 
astrophysics (improved satellite observations of the microwave background, analysis of su- 
pernovae and quasar radiation, space experiments for testing the WEP and constancy of the 
gravitational constant G , laboratary tests of FPC stability etc.) are expected to bring new 
results of major significance for the whole physical picture of the world. The discovery of 
a dynamical nature of a number of FPCs, as well as revision of the level of fundamentality 
for many of them, will inevitably lead to the corresponding revision in the foundations of 
theoretical metrology. Following the change in the general theoretical paradigm, the set of 
FPCs will also be revised: the parameters, whose dynamical nature will be proved, will lose 
their fundamental status, giving way to the basic constants of an underlying unified theory. 

If a does vary, a question of great interest for both physics and metrology is: which of 
the quantities in the definition (1) of a, namely, c, h, e oi Eq, are really varying. Different 
viewpoints may be found in the literature. 

The theories mentioned in Sec. 3 modify the Maxwell equations and actually introduce 
a scalar field dependence of the effective electronic charge e. On the other hand, Hehl and 
Obukhov [69] note that, since e and /i, as the units of charge and action, are invariants 
(4D scalars), to reconcile a variable a with the invariable (premetric) form of the Maxwell 
equations one may choose the speed of light c to vary. For variable c theories see, e.g., Ref. 
[70], the brief review [71] and references therein. Though, as noted by Alfonso- Fans [72], such 
first principles as local Lorentz invariance and local position invariance disfavour variations 
of c as compared to changes of e as a reason for varying a . 

Hehl and Obukhov [69] mention one more possibility: they represent a as 

a = e^no/iAnh) = l]o/(2i?i^), (27) 

where Qo is the impedance of free space and Rk is the von Klitzing constant (quantum Hall 
resistance). In these expressions, the speed of light c has disappeared. If, again, e and h 
are regarded real constants, the putative a variability should be attributed to the vacuum 
impedance Qq that becomes a dynamic (dilaton) field. 

Tobar [27] represented a in terms of the ratio of the quanta of electric ($e) and magnetic 
( $m ) fluxes of force of the electron: 

a = ^(-^] . (28) 
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The quantity a, being dimensionless, is independent of the system of units used. Eq. 
makes this independence manifest, which is an advantage as compared with the conventional 
expressions (1). The representation (j^H|) made it possible [27] to interpret the range vari- 
ation of a at the electroweak energy scale (see the end of Sec. 2.1) as appearing due to 
equal components of electric screening (polarization of vacuum) and magnetic anti-screening 
(magnetization of vacuum), which cause the perceived quanta of electric charge to increase 
at small distances, while the magnetic flux quanta decrease. Thus there emerge the concepts 
of "bare magnetic flux quanta" and "bare electric charge" . With regard to a putative drift 
at astrophysical or cosmological scales, Tobar [27] also interpreted it as a differential drift 
of the electric and magnetic fluxes, but concluded that it is impossible to determine which 
fundamental constant is varying if a varies. 

In our view, the constants c, related to the space-time structure, and h, related to the 
quantum properties of all kinds of matter, are more universal than e that characterizes the 
electromagnetic interaction only and are therefore more likely to be true constants. It is quite 
natural that unification theories modify the Maxwell equations by introducing slowly varying 
scalars which make the effective electronic charge vary. This view is also consistent with Hehl- 
Obukhov's and Tobar's interpretations. On the other hand, the constant Eq, which is present 
in the SI definition of the fine structure constant, only connects the electric SI units with the 
mechanical ones and does not have a fundamental meaning. (See also a recent discussion of 
SI and other systems of units in Refs. [73, 74].) 

Thus a variable a raises a number of questions in fundamental metrology, requiring 
further experimental and theoretical studies. 

As to the needs of practical metrology, the existing constraints on possible variations of 
a in the modern epoch, no more than a few units of the 17th significant digit per year [see 
Eq. (9)] , compared to the achieved measurement accuracy of 10^ , give no ground to expect any 
appreciable change in this constant in the foreseeable future. Thus there is no evident reason 
for taking into account the conjectural dynamical nature of a in the values of references 
of any physical quantities, their construction and analysis. However, as follows from our 
discussion, there are at least two exceptions: (i) observational and theoretical cosmology and 
astrophysics of remote quasars and galaxies, and (ii) high energy physics related to unification 
of interactions. 

A problem of importance for both theoretical and practical metrology is the connection 
between local and cosmological variations of the FPCs, a in particular. Such a connection 
should be known for comparison between, e.g., terrestrial and quasar measurements of a and 
a and their proper interpretation. On the other hand, since this connection appears to be 
model-dependent [52-54], spatial variations of FPCs deserve special attention due to their 
potential use for model selection in physics and cosmology. 
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